Compact stars in the QCD phase diagram II (CSQCD II) 

May 20-24, 2009, KIAA at Peking University, Beijing - P. R. China 

http : / / vega . bac . pku . edu . cn/ rxxu/ csqcd . htm 



Quark deconfinement phase transition in neutron 

stars 



Grigor Alaverdyan 
Department of Radiophysics 
Chair of Wave Processes Theory 
Yerevan State University 
Yerevan 0025 
Armenia 

Email : galaverdyanSysu . am 



Abstract 

The hadron-quark phase transition in the interior of compact stars is investi- 
gated, when the transition proceeds through a mixed phase. The hadronic phase 
is described in the framework of relativistic mean-field theory, when also the 
scalar-isovector delta-meson mean-field is taken into account. The changes of 
the parameters of phase transition caused by the presence of delta-meson field 
are explored. The results of calculation of structure of the mixed phase ( Glen- 
denning construction) are compared with the results of usual first-order phase 
transition (Maxwell construction). 
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1 Introduction 

The structure of compact stars functionally depends on the equation of state (EOS) 
of matter in a sufficiently wide range of densities - from 7.9 g/cm^ (the endpoint 
of thermonuclear burning) to one order of magnitude higher than nuclear saturation 
density. Therefore, the study of properties and composition of the matter constituents 
at extremely high density region is of a great interest in both nuclear and neutron star 
physics. The relativistic mean- field (RMF) theory [11 [21 [3] has been effectively applied 
to describe the structure of finite nuclei [H [5], the features of heavy- ion collisions [6i,i7|, 
and the equation of state (EOS) of nuclear matter [8]. Inclusion of the scalar-isovector 
(5-meson in this theoretical scheme and investigation of its infiuence on low density 
asymmetric nuclear matter was realized in Refs. [9l [TOl [11] . At sufficiently high den- 
sity, different exotic degrees of freedom, such as pion and kaon condensates, also 
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deconfined quarks, may appear in the strongly interacting matter. The modern con- 
cept of hadron-quark phase transition is based on the feature of that transition, that 
is the presence of two conserved quantities in this transition: baryon number and 
electric charge[T2]. It is known that, depending on the value of surface tension, cxs, 
the phase transition of nuclear matter into quark matter can occur in two scenarios 
|13j : ordinary first order phase transition with a density jump (Maxwell construc- 
tion), or formation of a mixed hadron-quark matter with a continuous variation of 
pressure and density (Glendenning construction) [12] . Uncertainty of the surface ten- 
sion values does not allow to determine the phase transition scenario, taking place 
in realty. In our recent paper [13] in the assumption that the transition to quark 
matter is a usual first-order phase transition, described by Maxwell construction, we 
have shown that the presence of the 5-meson field leads to the decrease of transition 
pressure Pq, of baryon number densities n^v and hq. In this article we investigate 
the hadron-quark phase transition of neutron star matter, when the transition pro- 
ceeds through a mixed phase. The results of calculation of structure of the mixed 
phase (Glendenning construction) are compared with the results of usual first-order 
phase transition (Maxwell construction). Influence of 5-meson field on phase transi- 
tion characteristics is discussed. Finally, using the EOS obtained, we calculate the 
integral and structural characteristics of Neutron stars with quark degrees of freedom. 



2 Neutron star matter equation of state 
2.1 Nuclear matter 

In this section we consider the EOS of matter in the region of nuclear and supranuclear 
density(n > 0.1 fm~^). For the lower density region, corresponding to the outer and 
inner crust of the star, we have used the EOS of Baym-Bethe-Pethick (BBP)[15j. 
For description of hadronic phase we use the relativistic Lagrangian density of many- 
particle system consisted of nucleons, p, n, and isoscalar-scalar (a), isoscalar- vector 
[uj), isovector-scalar (6), and isovector- vector (p) - exchanged mesons : 

Caujp5{(T{x),uj^{x),p^{x),6{x)) = C^ujf,{a{x),uj^{x),p^{x)) -U{cx{x)) + Cs{6{x)), (1) 

where Ccrujp is the linear part of relativistic Lagrangian density without (5-meson field 
[TBj . U{(j) = Im^r ((^o-O")^ + f (fi'o-O')'^ is the cr-meson self-interaction term and 

Cs{S) = gsi^NrNSiJN + lidjd^6- ms6') (2) 

is the contribution of the 5-meson field. This Lagrangian density (1) contains the 
meson-nucleon coupling constants, g„, g^, gp, gs and also parameters of a-field self- 
interacting terms, b and c. In our calculations we take for the 6 coupling constant. 
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0-5 = {gs/'msf' = 2.5 fm^, as in [TOl [TTl [T3], for the bare nucleon mass, rriM = 938.93 
MeV, for the nucleon effective mass, = 0.78 m^v, for the baryon number density 
at saturation, no = 0.153 fm"^, for the binding energy per baryon, /o = —16.3 
MeV, for the incompressibihty modulus, K = 300 MeV, and for the asymmetry 
energy, E^^l^ = 32.5 MeV. Five other constants = {g^/rriaY, '^'^ = {Quj / 'fn^jf , 
o-p = {9p/^p)'^j b and c then can be numerically determined [1^. In Table 1 we list the 





aa ( fm^) 


acu ( 


fm2) 


as ( fm^) 


ttp ( fm^) 


b ( fm-i) 


c 


aojp 


9.154 


4.^ 


^28 





4.794 


1.654- 10-2 


1.319 - 10-2 


aujp5 


9.154 


4.^ 


^28 


2.5 


13.621 


1.654- 10-2 


1.319 - 10-2 



Table 1: Model parameters without (crujp) and with {aujp6) a 6 -meson field. 



values of the model parameters with {aujp6) and without (crup) the isovector-scalar 
S meson interaction channel. 

The knowledge of the model parameters makes it possible to solve the set of 
four equations in a self-consistent way and to determine the re-denoted mean-fields, 

a = ga-a, ijj = QiJJo, S = gsS , and p = gp'po^^^ , depending on baryon number density 
n and asymmetry parameter a = {un — Up) / n. The standard QHD procedure allows 
to obtain expressions for energy density e{n, a) and pressure P{n, a) of nuclear npe 
plasma: 



k- (n,o) 
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£NM{n,a,pe 
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+ ^/Xe(/i^-mg'^'-^ J ^k^ + k'dk, (4) 



where He is the chemical potential of electrons, 

m*((T, 6) = rriN — cr — 5, 'm*^{a, 5) = rriN — a + 6 (5) 
are the effective masses of the proton and neutron, respectively, and 

k±{n,a)= i-^il±a)\ . (6) 
The chemical potentials of the proton and neutron are given by 

/ip(n, a) = \lkF {nf (1 - af^'^ + {rriN - cr - 5) ^ + u; + ^p, (7) 

fin{n,a) = y kp (n) (1 + a) ' + (mjv - + 5) +co'--p. (8) 




Figure 1: Energy per baryon func- 
tion of the baryon number density n and 
the asymmetry parameter a in case of a /3 
-equiUbrium charged npe-plasma. 
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Figure 2: Re-denoted meson mean-fields 
as a function of the baryon number den- 
sity n in case of a /3-equiUbrium charge- 
neutral npe-plasma with and without b- 
meson field. 



In Fig 1 we illustrate the 3D-plot of the energy per baryon, E\, (n, a) = eNhi/n, 
as a function of the baryon number density n and asymmetry parameter a in case 
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of a /3 -equilibrium charged npe-plasma[17]. The curves correspond to different fixed 
values of the charge per baryon, q = (n^ — ne)/n = (1 — a)/2 — n^/n. The thick 
one corresponds to /3-equilibrium charge-neutral npe-matter. The lower and upper 
surfaces correspond to the " oujp" and " oujpb" models, respectively. Clearly, including 
a b -meson field increases the energy per baryon, and this change is greater for larger 
values of the asymmetry parameter. For a fixed value of the specific charge, the 
asymmetry parameter falls off monotonically with the increase of density. 

In Fig 2 we plotted the effective mean-fields of exchanged mesons, a, cj, p and 8 
as a function of the baryon number density n for the charge-neutral /3-equilibrium 
npe-plasma. The solid and dashed hues correspond to the " aujpS''' and " aujp" models, 
respectively. From Fig 2 one can see that the inclusion of the scalar-isovector virtual 
(5(ao(980)) meson results in significant changes of p and b meson fields. This can 
result in changes of deconfinement phase transition parameters and, thus, alter the 
structural characteristics of neutron stars. 

The results of our analysis show that the scalar - isovector 5-meson field inclusion 
leads to the increase of the EOS stiffness of nuclear matter due to the splitting of 
proton and neutron effective masses, and also the increase of asymmetry energy (for 
details see Ref.[T7j). 



2.2 Quark matter 

To describe the quark phase an improved version of the MIT bag model[T8] is used, in 
which the interactions between m, (i, s quarks inside the bag are taken into account in 
the one-gluon exchange approximation fT9] . The quark phase consists of three quark 
flavors M, s and electrons in equilibrium with respect to week interactions. We 
choose rriu = 5 MeV, = 7 MeV and = 150 MeV for quark masses and = 0.5 
for the strong interaction constant. 



2.3 Deconfinement phase transition parameters 

There are two independent conserved charges in hadron-quark phase transition: bary- 
onic charge and electric charge. The constituents chemical potentials of the npe- 
plasma in /3-equilibrium are expressed through two potentials, pl^^^^ and pif^^\ 
according to conserved charges, as follows 

_ _ (NM) (NM) _ (NM) 

In this case, the pressure P, energy density e and baryon number density n, are 

r r . ^ (NM) J (nM) J-. , (NM) (NM)-, , (NM) (NM)-, 

functions of potentials, /i^ ' and pli ': FnmKpI \ Pel )? ^nmIpI \ Pel )^ 

(,,(NM) (NM). 
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The particle species chemical potentials for wrfse-plasma in /^-equilibrium are ex- 
pressed through the chemical potentials jj^l^^^^ and jJ^l^^^^ as follows 



1 f,,iQM) r, ,(QM)\ 



1 / (QM) , (QM)\ /-|^^ 
f^d = f^s = ^[fJ'l +fJ'el (10) 

_ _ (QM) 

In this case, the thermodynamic characteristics are functions of chemical potentials 

The mechanical and chemical equilibrium conditions (Gibbs conditions) for mixed 
phase are 

(QM) _ (NM) _ (QM) _ (NM) _ , x 

PQMifJ'b, fJ'el) = PNM^fJ'b, ^-el)- (12) 

The volume fraction of quark phase is 

X = Vqm/{Vqm + Vnm), (13) 

where Vqm and Vnm are volumes occupied by quark matter and nucleonic matter, 
respectively. 

We applied the global electrical neutrality condition for mixed quark-nucleonic 
matter, according to Glendenning 



;i - X) [npifib, fJ'el) - neiflel)] 



2 11 

- riuifJ'b, fJ'el) - - ridifJ'b, fJ'el) - - risifJ'b, fJ'el) - rieif^el) 



0. (14) 



+x 

The baryon number density in the mixed phase is determined as 

n= (I-X) [llpiHb, IJ'el) + riniflb, fJ'el)] 
+ 2 X [nuifJ-b, fJ-el) + ridif^b, fJ'el) + risifJ'b, f^el)] , (15) 

and the energy density is 

e = {1 -X) l^pil^b, l^el) + £^n(/ife, fJ'el)] 
+X [^uifJb, fJel) + ^difJb, fJel) + ^sifJb, fJel)] + ^eif^el)- (16) 

In case of x = 0, the chemical potentials /i^ and /i^ corresponding to the lower 
threshold of a mixed phase, are determined by solving Eqs. (12) and (14). This 
allows to determine the lower boundary parameters P/v, En and n^. Similarly, we 
determine the upper boundary values of mixed phase parameters, Pq, Sq and rig, for 
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Model 


(fm-3) 


(fm-3) 


Pn 
(MeV/fm^) 


Pq 

(MeV/fm^) 


En 

(MeV/fm^) 


(MeV/fm^) 


BGOaup 


0.0717 


1.083 


0.336 


327.747 


67.728 


1280.889 


BGOaujpS 


0.0771 


1.083 


0.434 


327.745 


72.793 


1280.884 


BlOOaup 


0.2596 


1.436 


18.025 


471.310 


253.814 


1870.769 


B100aujp6 


0.2409 


1.448 


16.911 


474.368 


235.029 


1889.336 



Table 2: The Mixed phase threshold parameters without [crujp) and with {auipS) a S 
-meson field for bag parameter values, i? = 60 MeV/fm^ and B = 100 MeV/fm^. 



X = 1. The system of Eqs. (12), (14), (15) and (16) makes it possible to determine 
EOS of mixed phase between this critical states. 

Note, that in the case of an ordinary first-order phase transition both nuclear 
and quark matter are assumed to be separately electrically neutral, and at some 
pressure Pq, corresponding to the coexistence of the two phases, their baryon chemical 
potentials are equal, i.e., 

P'NM {Pq) = P-QM (Po) ■ (17) 

Such phase transition scenario is known as phase transition with constant pressure 
(Maxwell construction). 

Table 2 represents the parameter sets of the mixed phase both with and without 
5-meson field. It is shown that the presence of 5-field alters threshold characteristics 
of the mixed phase. For B — 60 MeV/fm^ the lower threshold parameters, un, Sn, 
Pn are increased, meanwhile the upper ones ng, eg, Pq are slowly decreased. For 
B = 100 MeV/fm^ this behavior changes to opposite. 

In Fig 3 and Fig 4 we plot the EOS of star matter with deconfinement phase 
transition for two values of bag constant, B = 60 MeV/fm^ and B = 100 MeV/fm^, 
respectively. The dotted curves correspond to pure nucleonic and quark matters 
without any phase transition, while the solid lines correspond to two alternative phase 
transition scenarios. Open circles show the boundary points of the mixed phase. 

In Fig 5 we plot the particle species number densities as a function of baryon 
density n for Glendenning construction. Quarks appear at the critical density tin — 
0.241 fm~^. The hadronic matter completely disappears at uq — 1.448 fm~^, where 
the pure quark phase occurs. Fig 6 show the constituents number density as a function 
of baryon number density n for B = 100 MeV/fm^, when phase transition described 
according to Maxwell construction. Maxwell construction leads to the appearance of 
a discontinuity. In this case, the charge neutral nucleonic matter at baryon density 
rii = 0.475 fm~^ coexisted with the charge neutral quark matter at baryon density 
712 = 0.650 fm~^. Thus, the density range rii < n < ^2 is forbidden. In the Maxwell 
construction case the chemical potential of electrons, pe, has a jump at the coexistence 
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Figure 3: EOS of neutron star matter with Figure 4: Same as in Fig 3, but for B 
the deconfinement phase transition for a ^0*-* MeV/fm . 
bag constant -B = 60 MeV/fm^. For com- 
parison we plot both the Glendenning and 
Maxwell constructions. Open circles repre- 
sent the mixed phase boundaries. 



pressure Pq- Notice, that such discontinuity behaviour takes place only in usual first- 
order phase transition, i.e., in the Maxwell construction case. 

3 Properties of hybrid stars 

Using the EOS obtained in previous section, we calculate the integral and structure 
characteristics of Neutron stars with quark degrees of freedom. 

The hydrostatic equilibrium properties of spherical symmetric and isotropic com- 
pact stars in general relativity is described by the Tolman-Oppenheimer-Volkoff (TOV) 
equations PU]: 

G (P + £)(m-F47rr3p) 

r2 l-2Gm/r ' ^ ^ 

47rr2£, (19) 

where G is the gravitational constant, r is the distance from the center of star, P(r), 
sir) are the pressure and energy density at the radius r, respectively, and m(r) is the 



dP 
dr 
dm 
dr 
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Figure 5: Constituents number density 
versus baryon number density n for B = 
100 MeV/fm^ in case of Glendenning con- 
struction. Vertical dotted lines represent 
the mixed phase boundaries. 



Figure 6: Same as in Fig 5, but for 
Maxwell construction. Vertical dotted lines 
represent the density jump boundaries. 



mass inside a sphere of radius r. For integration of the TOY equations it is necessary 
to know the EOS of neutron star matter in a form e{P) . Using the neutron star matter 
EOS , obtained in previous section, we have integrated the Tolman-Oppenheimer 
Volkoff equations and obtained the gravitational mass M and the radius R of compact 
stars (with and without quark degrees of freedom) for the different values of central 
pressure Pc- 

Fig 7 and Fig 8 illustrates the M{R) dependence of neutron stars for two values of 
bag constant -B = 60 MeV/fm^ and B = 100 MeV/fm^, respectively. We can see, that 
the behaviour of mass-radius dependence significantly differs for two types of phase 
transitions. Fig 7 shows, that for B = 60 MeV/fm^ there are unstable region, where 
dM/dPc < 0, between two stable branches of compact stars, corresponding to config- 
urations with and without quark matter. In this case, there is a nonzero minimum 
value of the radius of quark phase core. Accretion of matter on a critical neutron star 
configuration will then result in a catastrophic rearrangement of the star, forming 
a star with a quark matter core. The range of mass values for stars, containing the 
mixed phase, is [0.085Mq; 1.853Mq] for B = QO MeV/fm^, and is [O.997M0; I.78OM0] 
for B = 100 MeV/fm^. In case of Maxwellian type phase transition the analogous 
range is [O.2I6M0; I.828M0] for B = 60 MeV/fm^. From Fig 8 one can observe that 
when B = lOOMeV/fm^ MeV/fm^, the star configurations with deconfined quark 
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Figure 7: The mass-radius relation of 
neutron star with different deconfinement 
phase transition scenarios for bag con- 
stant B = 60 MeV/fm^. Open circles 
and squares denote the critical configura- 
tions for Glendenning and Maxwellian type 
transitions, respectively. Solid circles and 
squares denote hybrid stars with minimal 
and maximal masses, respectively. . 
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Figure 8: Same as in Fig 7, but for B = 
100 MeV/fm3. 



matter is unstable. Thus, the stable neutron star maximum mass is 1.894M0. Our 
analysis show, that the pressure upper threshold value for mixed phase is larger, then 
the pressure, corresponding to maximum mass configuration. This means, that in the 
model, considered here, the mixed phase can exist in the center of compact stars, but 
no pure quark matter can exist. 



4 Conclusion 

In this paper we have studied the deconfinement phase transition of neutron star mat- 
ter, when the nuclear matter is described in the RMF theory with 5-meson effective 
fields. We show that the inclusion of scalar isovector (5-meson field terms leads to 
the stiff nuclear matter EOS. The presence of scalar isovector 5-meson field alters 
the threshold characteristics of the mixed phase. For B = 60 MeV/fm^, the lower 
threshold parameters, um, en, Pn are increased, meanwhile the upper ones, uq, eq, 
Pq are slowly decreased. For B = 100 MeV/fm^ this behavior changes to opposite. 
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For EOS used in this study, the central pressure of the maximum mass neutron 
stars is less than the mixed phase upper threshold Pq. Thus, the corresponding 
hybrid stars do not contain pure strange quark matter core. 



The author would like to thank G.S.Hajyan and Yu.L.Vartanyan for useful dis- 
cussions on issues related to the subject of this research. 
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